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ABSTRACT 

Synthetic spectra covering the wavelength range 900A to 3000A provide an 
accurate fit, established by a xl analysis, to a combined observed spectrum of 
RW Sextantis. Two separately calibrated distances to the system establish the 
synthetic spectrum comparison on an absolute flux basis but with two alternative 
scaling factors, requiring alternative values of M for final models. Based on 
comparisons for a range of M values, the observed spectrum does not follow the 
standard model. Rather than the exponent 0.25 in the expression for the radial 
temperature profile, a value close to 0.125 produces a synthetic spectrum with 
an accurate fit to the combined spectrum. A study of time-series FUSE spectra 
shows that a proposed warped or tilted disk is not supported by the data; an 
alternative proposal is that an observed non-axisymmetric wind results from an 
interaction with the mass transfer stream debris. 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS5-26555, and the NASA-CNES-CSA Far Ultraviolet Explorer, which is 
operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 
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Subject headings: Stars:Novae, Cataclysmic Variables, Stars:White Dwarfs, Stars: 
Individual: Constellation Name: RW Sextantis 



Introduction 



Cataclysmic variables (CVs) are semi-detached binary stars in which a late main se - 
quence star loses mass onto a white dwarf (WD) by Roche lobe overflow (jWarnerl Il995l ). 
In non-magnetic systems the mass transfer stream produces an accretion disk with mass 
transport inward and angular momentum transport outward, driven by viscous processes. 
The accretion disk may extend inward to the WD; the outer boundary extends to a tidal 
cutoff limit imposed by the secondary star in the steady state case. If the mass transfer 
rate is below a certain limit, the accretion disk is unstable and undergoes brightness cycles 
(outb ursts) , and if above the limit, the accretion disk is stable against outbursts (jOsaki 



19961 ). The latter objects (w 
(NL) systems. As sh own by 



rich have no recorde d outburst of any type) are called nova- like 
Cowley et al. ( 1977 ). Beuermann. Stasiewskifc Schwopel (1992 ) 



(hereaf t er BE U) and iPrinja et al. (b()03h . RW Sex is a NL system ( see IWarnerl Jl995h and 



Lasotal (120011 ) for a detailed discussion). RW Sex has a Hipparcos (jPerryman et al 



19971 ) 



parallax and a separately calibrated distance (BEU), of importance in constraining the mass 
transfer rat e M. BEU provid e estim ates of the component stellar masses and the orbital 
inclination; Greenstein fc Okd (119821 ) provide an estimate of M = 1.0xlO~ 8 M m yr _1 while 



Wadd (119881 ) quotes a value of M = 3.OxlO- 9 M yr" 1 from IPattersonl (11984 ). A FUSE 
spectrum of RW Sex is available as well as HST and IU E spectra. FU SE spectra are im- 
portant in constraining the WD T e g- as well as M since the peaks of the radiation curves fall 
in the FUSE spectral range. 

NL systems are of special interest because they are e xpected to have an accretion disk ra- 
dial temperature profile given by an analytic expression (IFrank. King fc Raindll992l . eq.5.41) 
(hereafter FKR, cf eq. (2) below) which defines the so-called standard model; the expression 



includes the m ass transfer rate M as an exp 



the period gap (IHowell. Nelson fc Rappaport 



icit v ariable. In NL systems which are above 



200 lj) the accre tion disk dominates the system 



spectrum (with the exception of BK Lyn (IZellem et al.ll2009l )). Fitting a synthetic spectrum 
based on the analytic model (a proxy of the accretion disk temperature profile) to an ob- 
served spectrum potentially determines M. T his physical parameter is of basic im portance 
since it controls the evolution of CV systems (IHowell. Nelson fc Rappaportll200ll ). But the 
analytic expression also is an explicit funct ion of the WD mass, M w( j, and the mass must 
be determined independently. IWadd (119881 ). using IUE spectra, showed that NL systems 
systematically disagree with the standard model when either black body spectra or Kurucz 
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stellar model spectra are used to represent the accretion disk. The Kurucz spectra were too 
"blue", i.e., had too large a spectral gradient as compared with the IUE spectra and this 
was specifically true of RW Sex. Wade concluded that the proble m was w i th use of Kurucz 
spectra; they are not suited to representation of accretion disks. iHubenyl (119901 ) has devel- 
oped a procedure for modeling annuli of accretion disks that explicitly includes calculation 
of synthetic spectra using the standard model. In this paper we apply the Hubeny model to 
determine system parameters for RW Sex and to test whether the system conforms to the 
standard model. 



2. The FUSE and HST/GHRS Spectra 



Table 1 prese nts the observations log. This is the same observation set studied by 



Prinja et al.l (120031 ) and discussed by those authors. 



2.1. The FUSE Spectrum 

The FUSE spectrum of RW Sex consists of 25 exposures (spacecraft orbits) totaling 
more than 25ks of calibrated (good) exposure time. The FUSE spectrum (flux and errors) 
were extracted using the CalFUSE software, initially binned at 0.016A. The data were re- 
binned at an interval of 0.5A using the IRAF command trebin. Typical noise levels are of 
order 3% of the signal level. Each spectrum consists of a list of flux values (erg cm~ 2 s _1 A ) 
and associated noise values in the same units. The FUSE spectrum covers the interval 
903A to 1188A. 

In Figure 1 we present the combined spectrum (25ks) annotated with absorption lines. 
It is impractical to place individual error bars on the plotted points. Besides the obvious 
Ly/3 and Ly7 broad features (due to the disk and accreting WD), the system exhibits broad 
absorption lines (which appear blue shifted in this composite by as much as 3A) and a 
multitude of very sharp absorption lines of interstellar origin. We describe these absorption 
lines below and list some of them in Table 2 and Table 3. The observed wavelengths of the 
sharp absorption lines listed in Table 2 were measured in the ISiC channel for A < 11 00 A , 
and 2LiF channel for A > HOOA ; as a consequence the observed ISM lines are blue shifted 
by ~ 0.05A for A < HOOA and red shifted by about the same amount for A > HOOA as the 
channels are not precisely aligned. When we co-added the channels and different orbits we 
aligned all the spectral segments. 

NL system spectra typically exhibit strong absorption lines from high excitation ions, a 
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signature of a wind believed to originate f rom the accretion disk ( IPereyra. Hillier fc Turnshek 
2006 ; Proga 2003al lbl). See the review by Drew fc Proga ( 2000 ). The presence of these lines 



is the largest current obstacle to spectrum synthesis modeling of the accretion disk in NL 
systems. In the present instance it is necessary to mask appreciable sections of the FUSE 
spectrum because we have no current ability to model the wind. This requirement reduces 
the number of lines that we fit in the FUSE spectrum. 



2.2. The ISM Hydrogen lines 

The FUSE spectrum of RW Sex exhibits a forest of molecular hydrogen absorption lines 
identified by the Werner (W) and Lyman (L) bands, upper vibrational level (ranging from 
16 to 1) and rotational transition R, P and Q with lower rotational state J=l,2 and 3. For 
example, starting at 933. 2A there is the L16R2 line, and next to it at 933. 6A the W4Q3 line. 
There are more than 100 molecular hydrogen lines extending to ~1100A (L0R2 1110. 10A). 
We have marked th ese lines (vertical tick mark) in the upper part of each panel in Figure 1 



(annoted MH). See ISembach et al.l (120011 ) for the entire list of molecular hydrogen lines and 
their wavelengths. 

The atomic hydrogen lines (Lyman series) are marked below each panel. 

The feature at 1152 A is a fixed pattern noise (FPN) due to the FUSE detectors. 



2.3. The ISM metal lines 

In addition to the hydrogen atomic and molecular lines, we identify many orders of 
neutral oxygen (Ol) starting at 916. 9A (the 26th order) and extending to about ~ 1000A 
with O I 4th and 3rd orders at AA1025.9 and AA1039.3. Not all the lines are listed in Table 2, 
but lines from all the orders (between 24u and 3u) are identified. In the comments column, 
lett er u stands for an ultr aviolet multiplet and the preceeding number is the multiplet number 



(see iMortonl (120001 . 120031 )). The oxygen lines are marked with vertical ticks (at mid- height; 
annoted OI). 

The FUSE spectrum of RW Sex is rich in oxygen lines in addition to metal lines which 
are frequently observed in FUSE spectra of CVs. These additional metal lines are An (AA 
1048.2, 1066.7), Cll (AA 1036.34, 1037.1), Fell (bottom panel of Figure 1), Nl (A -1034), 
Nil (A ~1084). We also identify many phosphorus (Pll) lines between ~ 960A and ~ 965 A 
as well as a single line at 1152. 8A (distinct from the well-known fixed pattern noise - FPN 
- introduced by the FUSE detectors at 1152A ; both this Pll line and the FPN feature are 



-5 - 



well documented (jSembach et al.l l200lh ) . It is not unusual to find that a FUSE spectrum 
rich in ISM 1 lines also exhibits many P II lines; P II lines are an excellent neutral oxygen 
tracer in various physical environments in the ISM and are particularly useful as a proxy for 
Ol lines when these are saturated or blended (jLeBouteiller et al.l 120081 ). 



2.4. The HST spectra 

The composite HST spectrum of RW Sex shown in Figure 2 is a combination of two 
GHRS spectra, z37v0107t and z37v0108t. 

The HST spectrum of RW Sex is characterized by deep and broad absorption lines from 
the CV source together with much shallower and sharper lines from the ISM. The lines are 
identified in Table 4. 

The broad absorption lines can be divided into two distinct groups. The first group of 
broad lines includes Sim (~ 1140, ~ 1205, ~ 1300, & 1327 A), Si IV (~ 1392 & 1401 A), 
Cm (~ 1173 A), and Hell (~ 1638 A) lines, all blue-shifted by ~1.8-2.5A with a width 
of about 5-8A. The second group of broad lines includes Nv (~ 1234 & 1237 A) and Civ 
(~ 1543 A) and is distinct from the first group in that the N and C lines are blue shifted 
by as much as 5A and 6A respectively. The lines also appear to be broader with a width 
of more than lOA. It is likely that the N v and C IV lines form in a hotter region expanding 
faster than the lines from the first group. The Lya line stands apart in that it is blue-shifted 
by less then lA. 

The sharp absorption lines all have a width of about ~ 0.5A or less and are blue-shifted 
by only 0.2-0.5A. The lines are listed in Table 4 and we tentatively identify, as in the FUSE 
spectrum, some P II phosporus lines at 1142.5, 1152.6 (also detected in the FUSE spectrum), 
and 1301.7 A (which could be affected by the Ol 1302 line). 



3. The IUE spectra 

Table 5 lists the low dispersion IU E spectra. There are three pairs of spectra covering 
the range 1150A to 3350A. We designate them as easel (SWP01671+LWR01583), case2 
(SWP02494+LWR03071), and case3 (SWP07500+LWR06494) in temporal sequence. Sep- 
arately, the SWP spectra cover the range 1150 A to 1978A and the LWR spectra cover the 
range 185lA to 3350A, tabulated in erg cm -2 s _1 A . The tabular interval of both spectra 
is about 2A and each has an associated tabulated noise level in the same units as the flux. As 
with the FUSE spectrum, typical noise levels are 3% of the signal. We used the IUEDAC 
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software package to extract the spectra, with zero correction for reddening, and including 
the a value at each t abular wavelength, and pro cessed the spectra with the standard re- 
calibration procedure (IMassa fc Fitzpatrickl |2000| ) . We examined the overlap region of each 
SWP and LWR pair and found that the easel pair fitted very well without scaling either 
spectrum. 



4. Merging FUSE, HST and IVE spectra 

Our model simulations concern fitting the observed continuum. The presence of broad 
absorption lines of high excitation species in the FVSE spectrum, which our model does not 
simulate, requires masking of parts of the FV SE spectrum and so reduces the utility of the 
FVSE spectrum in choosing among possible models. In addition, the short spectral range 
of the FVSE spectrum makes a stringent test of the spectral gradient fit more difficult. 
However, the peaks of the emission profiles of both the WD and the accretion disk fall in the 
FVSE range and the sensitivity to different WD T e g values and accretion disk M values is 
much greater than in HST or IV E spectra. In view of those features, we combine the FVSE 
and IV E spectra to obtain the best constraint on model parameters. The issue now is to 
choose the optimum IV E spectrum. We do not include the HST spectra in our simulations 
because of the restricted wavelength range of those spectra. 

Figure 3 shows the fit of the various spectra. Identifications of the spectra are in the 
figure legend. The mean of HST z37v0108t and z37v0109t (Table 1) is not plotted; it is 
similar to the mean of z37v0104t plus z37v0105t. In this overlap region the IVE easel 
(orange line) agrees well with the HST (blue) spectra. The other two IVE spectra are 
roughly accordant, with the exception of the case3 La geocoronal A 1220 line. The case-to- 
case variation of the IVE spectra indicates some temporal variation that is most marked in 
the absorption lines. 

The black line is the FVSE spectrum as observed; it is discrepant from both the HST 
and IVE spectra. Two possible explanations of the discrepancy are: (1) there was an error 
in the reduction of the spectrum, and, (2) RW Sex is slightly variable and was at a lower 
luminosity level at the time of the FVSE observation. We have checked t he reduction 



procedure and find no error. The MAST preview and the iPrinja et al.l (120031 ) paper agree 
with our data tabulation. The mutual agreement of the HST and IVE spectra support the 
postulate that the FVSE spectrum was anomalously weak at the time of observation. The 
red line is the FVSE spectrum divided by 0.82; the fit to the HST spectra now is excellent 
and, within the noise level of the IVE spectra, the fit to the IVE spectra is good. 
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We deleted the region of IUE easel that overlaps the FUSE spectrum and added the 
remaining part of easel to the FUSE spectrum divided by 0.82. This combined spectrum 
shows strong absorption lines associated with a wind/chromosphere (Figure 1) and other 
features like the P Cygni CIV line at 1540A (Figure 2). The synthetic spectrum does not 
model those features and it is necessary to mask them. At the same time, the widths of 
strong absorption lines are discriminants among different models. We were careful to mask 
features that clearly arose from a wind/chromosphere while minimizing the total amount of 
masking. 



The FUSE data we simulate consists of flux values, Fi (erg cm 2 s 1 A 1 ), and asso- 
ciated noise values, a x . Calculation of the reduced xl follows from the equation 

v JW(V ' 



Xl = l/(N -M)J2 



where N is the number of observed wavelengths, M is the number of model parameters 
(identified and discussed below), Fi is the observed flux at a particular wavelength, Oi is 
the associated uncertainty of Fi, and /(A,) is the flux calculated by the model. Our model 
consists of calculated flux values at a tabular interval of O.lA; the fitting process interpolates 
among the calculated flux values to the exact tabular wavelength of each observed data point. 



A rough test of a good model fit to an observed data set is a Xu 1 (iTaylorl 11982 
ch.12.3). In our case, a final xl^ would indicate essentially a perfect fit at the noise level- 
an accuracy unachievable in the presence of unmodeled spectral features. However, if we do 
use the available noise level tabulation, the calculated xl variation from model to model still 
provides a measure of relative model-to-model quality of fit. The subsequent discussion, §8., 
includes a more detailed test. 



5. Initially adopted system parameters 



By the basic paradigm of CV evolution ( jHowell. Nelson fc Rappaportl 12001 ) there is a 
close relation between orbital period and donor mass. In a recent paper, iKniggd (120071 ) has 
calibrated the relations among P, M 2 , R2, T e ff t2 , and donor spectral type. Thus, given an 
orbital period below the tabular upper limit of 6 h , M2 can be determined independently of 
other parameters. If there is an observationally-determined mass ratio, the WD mass follows 
directly. 

Table 6 (and tablenotes) lists initially adopted parameter s and their sources. Several of 
the system parameters are poorly known. IKniggd (120061 120071 ) determines M2 = 0.67M Q for 
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P = 0.24507 d . BEU determine a mass rat io 1/q — 1.35±0.1 which, with M%, produces M w d = 
0.9. IPanei. Althaus fc Benvermtol (boooh lists a WD radius of 8.82xl(r 3 i? for a 0.90M© 
homogeneous zero temperature Hamada-Salpeter carbon model WD; in our subsequent study 
of the observed spectra we correct the radius for the adopted WD T eS . 



Verbuntl (119871 ) determined a preferred value of E(B — V) = 0.0 with a ro ugh upper limit 
of E(B - V) of 0.03, based on a study of the 2175A "bump" in I UE sp ectra, bruch fc Eneel 
f|l994h list E{B - V) = 0.02. iMauche. Raymond fc Cordova! fjl988h determine a value of 
E(B - V) = 0.014 with appreciable uncertainty. We adopt E(B - V) = 0.0. 

A number o f studies have considered the t idal cutoff boundary, r,j, of accretion disks 

( Paczynski ( 1977 ); Papaloizou &: Pringle ( 1977)) ; Whitehurst ( 1988 ); Schwarzenberg-Czerny fc Rozyczk 
( Il988[ ): Iwhitehurst fc KineJ ( Il99l[ ): iGoodmarJ (|l993j)). These authors agree on r d ~0.33L>, 



:a 



where D is the separation of the stellar components. We adopt this expression for the tidal 
cutoff radius of the accretion disk. 



The Hipparcos (jPerryman et al.l 119971 ) parallax, 3.46±2.44mas, corresponds to a dis- 
tance of 289pc. BEU determine a distance of 150pc, based on an application of Bailey's 
method (jBaileylll98ll ) and different from the Hipparcos value by slightly more than la. We 
will study our results in the context of both distance determinations and we initially adopt 
the Hipparcos value. 

With respect to the Table 6 parameters, the orbital period is known with essentially 
perfect accuracy as compared to the other parameters. Of the remaining parameters Mi 
has no effect on the model spectrum; however T e s of the WD does affect the model and we 
include it as an adjustable parameter. 



6. The analysis program: BINSYN 



Our analysis uses the progr am suite BINSYN (ILinnell &: Hubenylll996l ); recent papers 
(ILinnell et al.ll2007l . l2008ai 120091 ) describe its application to CV systems in detail. Briefly, an 
initial calculation produces a set of annulus models for a giv en WD mass, radius, and mass 
trans fer rate. This calculation uses the program TLUSTY[J ( Hubenyill988l ; iHubeny fc L 



anz 



19951 ). TLUSTY considers the radial and vertical structure of the disk independently; the 
radial structure is based on the standard model (FKR) and so follows the prescribed relation 
between local T P ff and th e annu lus r adius. The vertical s tructu re is solved, self-consistently, 
as described by IHubeny! (Il990i ) and IHubeny &: Hubenyl (119981 ). The set of annulus models 



http: / /nova. astro. umd.edu 
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covers the accretion disk from its innermost (WD) radius to r/r w d.o 
the radius of the zero-temperature WD. 



50.0, where r w do is 



The primary sou rce of viscosity in CV ac c retion disks is b elieved to be MRI (mag n etoro - 
tational instability) faalbus fc Hawlevi Il99ll iBalbusI liooih . iHirose. Krolik fc Stone! f|2006h 
calculate MHD models with local dissipation of turb ulence and show th at the vertical extent 
of an annulus is greater than in previous models. iBlaes et al.l (120061 ) show that magnetic 
support has a significant effect on synthetic spectra of black hole (BH) accretion disk annuli 
in the X-ray region and illustrate the effect in the case of a BH system with Mbh = 6.62M 
and with an adopted a = 0.02. There is an insigificant effect in the visible and UV. They 
are able to simulate the effect of magnetic support within the TLUSTY framework by ad- 
justing the TLUSTY parameters Co an d Ci an d use this simulation to calculate the effect 
on synthetic spectra. TLUSTY allows for a vertical viscosity profile within a given annu- 
lus; it introduces an assignable division point within an annulus, dividing the annulus into 
deep layers and outer layers. Within the deep layers, viscosity follows a power law varia- 
tion specified by Cn , while the outer layers follow a power la w variation specifi ed by (j. See 
Hubeny fc Hubenyl (119981 ) for details. The primary effect (IBlaes et al.l I2006T) of magnetic 
support is a hardening of the spectrum in the 0.45keV region. iKing. Pringle fc Livid (120071 ) 
point out that MHD models that produce viscosity via MRI require a values that are a 
factor of 10 smaller than the a=0. 1-0.4 required by observational evidence and suggest that 
some caution still is needed in accepting the MHD results. Our models have used the default 
values Co = Ci = 0-0; this choice produces a constant viscosity vertically within the annulus. 
We have calculated a synthetic spectrum of one annulus that contributes significantly to the 
system synthetic spectrum using the same viscosity parameters used by IBlaes et al.l (120061 ) 
and have compared that synthetic spectrum with our model synthetic spectrum for the same 
annu lus for the spectra l range 900A to 3000A. We find no detectable difference, in agreement 
with lBlaes et al.l (120061 ). Magnetic supp ort is believed to extend an accretion disk verti cally 



as compared with the standard mode l (IHirose. Kro 



effect confirmed in the case of IX Vel (ILinnell et al 



ik fc Stone! 120061 : IBlaes et all 12009 ). an 



20071 ) . The absence of an eclipse in RW 



Sex prevents a test in this system. We assume a fixed a for the entire accretion disk. 

Table 7 lists properties of annuli calculated with TLUSTY(v.203n) for a mass transfer 
rate of 5. Ox lO~ 9 M /yr _1 . This illustrative case is from among the cases we tested in our 
simulations. The TLUSTY control file to calculate a given annulus requires a radius of 
the WD in units of Rq. All of the annuli used the radius of a zero temperature WD for 
a hom ogeneous carbon Hamada-Salpeter O.9OM model from iPanei. Althaus fc Benvenuto 
(120001 ). All of the annuli are solar composition models, and the models through r/r W d,o = 
26.00 are converged non-LTE models. The remaining models are so-called grey models (see 
the TLUSTY Users Guide for an explanation). The annulus spectral flux levels vary by 
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a factor of order 4 dex from the innermost annulus to the rim annulus. The grey models 
contribute very little to the accretion disk synthetic spectrum and the difference between the 
contribution of a grey model synthetic spectrum and a corresponding converged non-LTE 
model spectrum is negligible. Our adoption of solar composition models implicitly assumes 
that the secondary star, which supplies the material for the accretion disk, has an atmosphere 
with a solar composition. 

The annulus calculations adopted C,Mg,Al,Si, and Fe as explicit ions in addition to H 
and He. Tests show that there are detectable differences between synthetic spectra for only 
H,He as explicit ions and models including the metals listed. The remaining grey model 
annuli use the same set of explicit ions (see the TLUSTY Users Guide for an explanation). 
We believe this is the first instance in which non-LTE annulus models using metals as 
explicit ions have been used in a CV simulation. As we discuss subsequently, this study 
performed a xl analysis with M values of 1.0, 2.0, 4.0, 5.0, 5.25, 5.5, 5.75, 6.0, 6.25, 6.5 and 
7.0 xlO -9 M Q yr _1 . Each value of M has an associated table similar to Table 7, for the same 
range of annulus radii, with non-LTE models through r/r w d,o — 26.0 (the limit of the range 
of convergence) and for the explicit ions listed above. 



A IShakura fe Sunyaevl (119731 ) viscosity parameter a = 0.1 was used in calculating all 
annuli. Each line in Table 7 represents a separate annulus. Temperatures are in Kelvins. 
The column headed by mo is the column mass, in g cm -2 , above the central plane. The 
columns headed by z and N e are, respectively, the height (cm) above the central plane for 
a Rosseland optical depth of ~0.7 and the electron density (cm -3 ) at the same level. The 
column headed by log g is the log gravity (cgs units) in the z direction at a Rosseland optical 
depth of ~0.7. The column headed by tr oss is the Rosseland optical depth at the central 
plane. We call attention to the fact that the annuli are optically thick to the outer radius of 
the accretion disk (and true for the full range of M values). 

Following the calculation of annulus mo dels for assigned M w d and M, program SYN- 
SPEC(v.48) □ (IHubeny. Lanz fe Jefferylll994j ) was used to produce a synthetic spectrum for 



each annulus, at a spectral resolution of O.lA. We adopted solar composition for all synthetic 
spectra. The synthetic spectra include contributions from the explicit ions listed above and 
the remaining first 30 atomic species. 

The xl analysis mentioned above needs WD synthetic spectra to pair with each M ac- 
cretion disk model. We calculated non-LTE WD synthetic spectra for T e g- values of 35,000K, 
40,000K, 45,000K, 50,000K, 55,000K, 60,000K, 65,000K, 70,000K, and 75,000K with the 
same list of explicit ions listed above and at a spectral resolution of O.lA. This WD repre- 



2 http://nova. astro. umd.edu 
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sentation assumes that the solar composition mass transfer stream, after passing through 
the accretion disk to deposition on the WD, so adulterates the WD photosphere that it 
preserves a solar composition in spite of the inward diffusion of heavy ions through the WD 
atmosphere. 

BINSYN models the complete CV system, including the WD, secondary star, accretion 
disk face, and accretion disk rim as separate entities. The model represents phase-dependent 
and inclination-dependent effects, including eclipse effects and irradiation effects, on all of the 
system objects. Calculation of synthetic spectra for the stars requires polar T eS values and 
gravity-darkening exponents. We adopted a standard gravity-darkening exponent of 0.25 for 
the WD; the secondary star makes a negligible contribution to the system synthetic spectrum, 
as shown by a comparison of the secondary star synthetic spectrum and the accretion disk 
synthetic spectrum, and we do not list the secondary star simulation parameters. BINSYN 
represents the accretion disk by a specified number of annuli (45 in the present case), where 
that number typically is larger than the number of TLUSTY annulus models (22 in the 
present case; see Table 7). BINSYN calculates a synthetic spectrum for each of the annuli 
specified in the BINSYN system model by interpolation among the array of TLUSTY annulus 
models, with proper allowance for orbital inclination. The accretion disk T e g profile may 
follow the standard model, but, alternatively, the profile may be specified by a separate 
input file or by a non-standard radial temperature gradient as discussed below. 

BINSYN has provision to represent a bright spot on the rim face, but our present data 
do not require this facility (§9). 

7. Standard model simulations of RW Sex 

We adopted system parameters listed in Table 6 for an initial test; the initial value of M 
was arbitrary and was meant to test whether its value produced a synthetic system spectrum 
with even an approximate fit to the observed spectrum. On the assumption that the WD 
contribution would be small, we adopted a WD T eff = 50, 000K based on our previous studies 
of CV systems. If necessary, new models could be calculated with a revised WD T e fj. The 
initial choice of M (M = 1.0 xlO -9 M & yr -1 ) was obviously too small; we successively tried 
M = 2.0, 4.0, and 5.0 xlO -9 M & yr _1 and the last choice indicated a value that might be 
appropriate but with a problem described below. 



In constructing a model accretion disk it is necessary to allow for the change in WD 
radius from the zero temperature model. We used Table 4a of iPanei. Althaus fc Benvenuto 
feoooh to estimate the 50,000K WD radius at O.OO99i? . Our BINSYN model specified 45 



- 12 - 



annuli for the accretion disk, with assigned radii and corresponding standard model T e g 
values listed in Table 8. The first column in Table 7 and Table 8 specifies inner radii of 
annuli measured in units of the WD radius; in Table 7 the unit is a zero temperature WD 
and in Table 8 the unit is the WD radius for its assigned T eS of 50,000K, a larger quantity. 
Thus in Table 8 a given number in column 1 corresponds to a larger physical distance, and a 
lower annulus T e g, than for Table 7. The synthetic spectrum for a given BINSYN annulus is 
calculated by interpolation (temperature-wise) among the Table 7 entries. Ideally, we could 
recalculate the whole series of annulus models (via TLUSTY) based on the new value of 
R w d, but the change in the individual annulus models would be very small and would not 
be warranted in view of other uncertainties (e.g., the system parallax). 

Table 7 ends at a T e g of 6970K. The net effect is that the calculated BINSYN outer 
annuli have synthetic spectra corresponding to 6970K even if the specified standard model 
calls for a lower temperature. The effect on the synthetic spectrum of failing to follow 
the standard model is negligible. Adding a lower temperature (3500K) stellar synthetic 
spectrum to the array of Table 7, thereby permitting interpolation to a lower temperature in 
the outer BINSYN annuli, made no detectable change in the system synthetic spectrum. As 
separate issues, the outer annuli temperatures should not fall below 6000K to avoid outbursts 



(IOsakilll996l ; lLasotall200ll ) and tidal dissipatio n is expecte d to raise the outer accretion disk 



temperature above the standard model value (jSmak 



2002( 1. 



The adopted Hipparcos distance provided a fixed scaling factor by which to divide the 
synthetic specrum to superpose it on the observed ("corrected" FUSE plus IUE easel) 
spectrum. The spectral gradient of the superposed synthetic spectrum was much too large 
but the total integrated flux, judged by an eye estimate of the area below the synthetic 
spectrum, appeared about right. We repeated the entire process with M values, including 
the initial test, of 5.0, 5.25. 5.5, 5.75, 6.0, 6.25, 6.5, and 7.0 xl0~ 9 M yr" 1 . In all cases 
we used a WD T e g of 50,000K; the spectral gradients of all models were too large. Figure 4 
shows the results for five models, identified in the figure legend, and Figure 5 shows the 
results in the FUSE spectral region. The usable range of the IUE spectra extend to 3000A, 
matching the spectral interval covered by our synthetic spectra. 

Smaller mass transfer rates would produce an accretion disk with too small luminosity 
to fit the observations (most of the radiation curve would lie below the observed spectrum) 
and larger mass transfer rates would produce an accretion disk with too large luminosity, 
based on the adopted system distance. The contribution of the WD cannot be the source of 
the too large calculated spectral gradient because its contribution is much too small; a plot 
of the accretion disk contribution alone shows nearly identical results. On a purely empirical 
basis we inquire whether an analytic model with a smaller than standard radial temperature 
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gradient but the same total luminosity could fit the observations. 



The following equation defines the standard model (IFrank. King fc Raindll992f ). 



(3GMM 



a 



R* \ 1 I 2 



R 



1/4 



where M is the mass of the WD, M is the mass transfer rate, and R* is the radius of the 
WD. G is the gravitation constant and a is the Stefan-Boltzmann constant. A smaller radial 
temperature gradient but the same flux follows from 
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T eS (R) = SCL\ l-(^) ) , (3) 



R 



X 



EXP 



where EXP is a number smaller than 0.25 and SCL is an empirically-determined number 
to preserve the standard model T e s on the left side of the equation (and thereby preserve 
the same total flux since flux = o"T c 4 ff ). The significance of this generalization is that EXP 
becomes an additional model parameter (M in equation 1 now becomes 6). 

Equation 3 was included in BINSYN at the time of its development but this is the first 
instance of its use in a simulation of a CV system. See §10 for further discussion. 

The calculation proceeds as follows; an initial run at EXP = 0.25 (SCL = 1.0) deter- 
mines the standard model accretion disk luminosity which follows from Lum = J^li vT^^A 



where a is the Stefan-Boltzmann constant, T e fj; is the effective temperature of the j th an- 



j is the effective temperature of the j th 
nulus, Aj is the radiating area of the annulus, and there are M annuli. Then, having chosen 
a different value of EXP, repeated runs with modified values of SCL eventually identifies 
the SCL value producing the same accretion disk luminosity as for EXP = 0.25. During 
each run, application of equation 3 establishes the T e g values at the various annulus radii as- 
signed by BINSYN. Interpolation among the TLUSTY annulus spectra (Table 7) produces 
the BINSYN annulus spectra and integration over them produces the accretion disk syn- 
thetic spectrum. Since the temperature gradient is smaller for assigned EXP values than 
the standard model it is certain that the TLUSTY models will span the required BINSYN 
range of models. We found that the quality of fit improved as EXP is reduced until we 
reached the value EXP = 0.125 but then became poorer as EXP approached 0. Figure 6 
shows the results over the full spectral interval; Figure 7 shows the results in the FUSE 
region. 

To study the quality of fit in detail we have performed a xl analysis with 8 values 
of M (5.0, 5.25, 5.5, 5.75, 6.0, 6.25, 6.5, and 7.0 xl0~ 9 M & yr" 1 ) and 11 values of EXP 
(0.25, 0.225, 0.20, 0.175, 0.15, 0.125, 0.10, 0.075, 0.05, 0.025, and 0.0). As explained in 
connection with the <Tj of equation 1, our use of noise level <7j errors would guarantee that 
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the calculated xl values will be far larger than 1.0; nevertheless the calculation preserves 
relative values enabling us to identify best fit parameters. In the overlap region between the 
FUSE spectrum and the IUE spectrum the FUSE a\, of order lxl0~ 13 erg cm -2 s _1 A , 
are approximately a factor 2 smaller than the IUE a\. To give equal weight to all the 
observed points we used the mean a x of both data sets as a fixed o\ in equation 1. The 
calculated xl values, each with 1086 degrees of freedom, are in Table 9 and a plot is in 
Figure 8; standard model accretion disks are represented by the horizontal line at the top 
of the plot. For ease of calculation we have identified the degrees of freedom with N in 
equation 1; the effect is a very slight reduction of xl as compared with use of N-M as the 
degrees of freedom; it is equivalent to masking an additional 2A in the FUSE spectrum. A 
separate calculation for M = 5.5xlO _9 M yr _1 and EXP = 0.125 and using the individual 
FUSE and IUE o\ values in equation 1 produced a xl value of 18.7 in contrast to the value 
listed in Table 9. Use of the individual a\ in calculating Table 9 would have produced a 
Figure 8 essentially identical to the one plotted but with slightly larger contour labels. We 
call attention to the fact that the best fit value of EXP is close to 0.125 at all mass transfer 
rates, with a slight downward slope with increasing M. It is curious that the value 0.125 is 
half of the standard model value, 0.25. 

We used a WD T c g- = 50, 000K for all models appearing in Table 9. Although the 
WD contribution to a system synthetic spectrum is small, it is of interest to verify that the 
adopted T e s produces optimum xl fits. Further, having identified best fit values of EXP 
across the full range of M values tested, it is of interest to determine whether further refine- 
ment can identify a preferred value of WD T e g. We adopt the value EXP = 0.125 as fixed 
and repeat a xl analysis with the same values of M as Figure 8 and with 9 values of T e g 
35,000, 40,000, 45,000, 50,000, 55,000, 60,000, 65,000, 70,000 a nd 75,000K. Since the WD ra- 



dius is temperature-sensitive we determined 0.90M Q WD radii (IPanei. Althaus fc Benvenuto 



20001 ) for the various assigned T e g values and used them in calculating the accretion disks. 
The xl values are in Table 10 and a plot of the results is in Figure 9. It is clear that 
the choice of T e g = 50, 000K was about optimum for the available data but that no actual 
determination of the WD T e g is possible. 



8. The orbital inclination, distance, and a best fit model 

The Table 9 and Table 10 results apply for the adopted values of distance (289pc) and 
i (34°). 

According to BEU the value of i is 34°±6°. We have tested the effect on the synthetic 
spectrum fit, for M = 5.75xl0 -9 M & yr _1 , EXP = 0.125, by calculating a synthetic model 
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spectrum for the limiting values of i with all other parameters fixed. For i = 28°, as expected 
from the angle cosine, the synthetic spectrum was higher by 7%, and for i = 40° the synthetic 
spectrum was lower by 7%. These limiting spectra differ from the observed spectrum by a 
small amount, far smaller than the variation from the uncertainty in the distance, and could 
be compensated by a small change in the adopted distance without affecting the quality of 
fit. Our results are consistent with % = 34° but the synthetic spectra cannot constrain i. 

The Hipparcos parallax has la limits of ±2.44mas; the corresponding upper and lower 
distances are 980pc and 169pc. The RW Sex distance according to BEU is 150pc. If the 
true distance is close to 150pc a smaller synthetic spectrum scaling factor applies and a 
correspondingly smaller M is required. Experiments with smaller M values show that the 
same standard model spectral gradient problem persists in those cases as with the simulations 
discussed above. After several tests the value M = 2.0 xlO -9 M & yr _1 appeared the best fit 
using the 150pc scaling factor. Models with a range of values of EXP, equation 3, led to a 
visually identified best fit with EXP = 0.15, SCL = 0.51443. It is least laborious to make 
final adjustments in the scaling factor rather than recalculate models for slightly modified 
M; our final scaling factor is 2.5 xlO 41 corresponding to a distance of 162pc and in agreement 
with the determination of 150pc. The corresponding xli f° r comparison with Table 9, was 
15.61, indistinguishable from the optimum value from Table 9. For comparison, the scaling 
factor for Figures 4,5,6, and 7 is 9.14xl0 42 . 

We call attention to the fact that the smaller spectral gradient associated with EXP = 
0.15, as compared with the standard model, avoids the problem of disk instability against 
ouburst near a T e g of 6000K. The rim temperature of our best fit model is 7667K while the 
rim temperature of a corresponding standard model would be 4216K, well below the limit 
for onset of outbursts. 



Using the K-band 2MASS ( Skrutskie et a 
Mk magnitude for RW Sex from 



. 2006) magnitude of 10.07 and the absolute 



Knigge! J2006I . l2007h for an orbital period of 5.88 h , we derive 



a lower limit distance of 12 4pc and an upper limit of 217pc, slightly favoring the distance of 



150pc. ICowley et al.l (119771 ) quote a trigometric parallax determination of 7mas by Osvalds 
corresponding to a distance of ~143pc. The xl values do not distinguish between the 289pc 
and 150pc distances; the weight of the evidence favors 150pc which we accept. 

Our best fit model has the parameters listed in Table 11. 

Figure 10 shows the best fit model together with the various data sources, identified in 
the figure legend. Figure 11 shows the same data set for a restricted wavelength range. The 
observational data plotted are the values before masking, in contrast to all previous plots. 

Examination of Figure 11 shows that the model fits Lyman a fairly well although the 
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model has a wider and less deep profile than the HST spectrum. The FUSE Lyman f3 
line shows emission peaks in both wings relative to the model. We suggest this proposed 
unmodeled effect (emission wings) may be the explanation of those residuals and that the 
same explanation may apply to the Lyman 7 line. An extension of the argument could 
explain the remaining shorter wavelength residuals. 

Figure 12 is a plot of the residuals from the best fit model to the masked combined 
spectrum for a restricted wavelength range. The plot is a continuous line but the plot data 
consist of individual values at spacings of approximately 1.5A. Comparison with Figure 10 
shows that the synthetic spectrum is much smoother; Doppler shifts in the Keplerian accre- 
tion disk wash out narrow line features in the model spectrum; the residuals in the upper 
panel of Figure 10 have the appearance of noise, with a visually estimated mean of about 
2xl0 _13 erg cm -2 s _1 A and frequent larger residuals. The bottom panel of Figure 12 
shows that the residuals have a systematic variation in the 900A to 1200A interval. Ex- 
cluding that interval temporarily, which includes the entire FUSE spectrum, we test the 
hypothesis that the model spectrum is a good fit to the (masked, IUE) observed spectrum. 
There are 826 data lines in the spectral interval selected; the standard deviation, a T , of the 
residuals is 2.22 xlO -13 erg cm -2 s _1 A , the value of x 2 1S 828.62 and the value of xl is 
1.009 for 821 degrees o f freedom. Using the incomplete gamma function routine GAMMQ 



from iPress et al.l (119861 ) we calculate the pro bability Q t h at the x value as poor as 828.62 



could occur by chance is 0.419; according to IPress et al.l ( 119861 ) a value Q greater than 0.1 



indicates an acceptable model. We also apply the Kolmogorov-Smirnov test where the null 
hypotheis is that the residuals represent a normal distribution with mean of zero. The test 
determined the a r listed above and sorted the residuals into 60 bins, 58 of width 0.1cr r , one 
containing all residuals greater than 2.55cr r and the other containing all residuals smaller 
than — 2.55<r r . The program calculated a theoretical cumulative distribution function, the 
actual cumulative distribution function of the residuals, and the difference between the two. 
The maximu m difference was 0.168 while the critical value at the 0.05 level of significance 



( lOstld Il969l . Table 16) was 0.0473; the hypothesis (normal distribution with mean of zero) 



is rejected. 

The maximum systematic residuals of the FUSE spectrum, Figure 12, bottom panel, 
are of order 1/3 of the observed spectrum flux values. The best fit model, essentially identical 
with the magenta line of Figure 7, falls below the FUSE spectrum shortward of 1050A and 
has greater flux longward. We have already speculated that there may be unmodeled emission 
in the shorter wavelengths. We found (§4.) that the FUSE spectrum apparently was 
anomalously weak at the time of observation and a correction factor was necessary to combine 
it with the IUE spectrum. 
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8.1. How robust is the best fit model? 

The model parameters subject to adjustment are M w a, T e fj jW d, i, M, distance, and EXP. 
As explained above, variation of i produces a minor effect at the adopted i and in any case 
i cannot be used to constrain the model because RW Sex is a non-eclipsing system. 

Variation of M W( j has two effects; first, there is a variation of WD radius with mass, 
leading to variation of the WD luminous flux for a given WD T e ff, and second, the variation 
of the depth of the potential well produces variation of annulus T c g at a given distance 
from the WD for a given M. As shown by, e.g., Figure 10, the WD contribution to the 
system synthetic spectrum is small and the variation of it (the synthetic spectrum) due to 
a change in M w( j is second order; we neglect it in comparison with other major contributors 
to the model spectrum. Variation of an annulus T e g from variation of M w d can be included 
in the study of M effects since the latter also produce variation of annulus T e g at a given 
distance from the WD and the variation with M is the dominant effect. We exclude M w< j as 
a significant contributing source to the discrepancy between observation and the standard 
model. 

If we reset EXP = 0.25 (i.e., assert the standard model is applicable), can we vary 
the RW Sex distance or M or the WD T eff and achieve as good a fit to the observed spec- 
tra as our model does in Figure 10 and Figure 11? From Figure 4, the closest synthetic 
spectrum fit to the observations is with the smallest M, 5.0xlO _9 M Q yr _1 . Still smaller 
M values, tested through l.OxlO _9 M yr _1 , show the same trend and ultimately the M 
value becomes small enough that the accretion disk becomes unstable against outburst. At 
M = 2.OxlO _9 M yr _1 the outer accretion disk rim has a calculated T e g of 4215K, making it 
susceptable to outbursts. However, stream impact and tidal effects heat the outer disk region 



(jLasotal 120011 ) ; let us assume these effects stabilize the standard model disk (and ignore the 
consequent departure from the standard model) and inquire how well the standard model for 
M = 2.OxlO" 9 M yr" 1 fits the observations. Fi gure 13 shows the fit of the standard model 
to the observations. The system synthetic spectrum, including the contributions of both the 
accretion disk and the WD, has been divided by 3.0 xlO 41 , corresponding to a distance of 
180pc. A larger divisor shifts the synthetic spectrum downward, producing larger long wave- 
length discrepancies, and a smaller divisor shifts the synthetic spectrum upward, producing 
larger short wavelength discrepancies. The contribution of the 50,000K WD is at the bottom. 
Removing its contribution entirely would only marginally reduce the discrepancy from the 
system synthetic spectrum. A proposed higher T e s WD would raise the question of its source 
of heat but, in any case, the rapidly increasing WD FUV flux would only make the system 
synthetic spectrum discrepancy worse. Still smaller M values exhibit the same discrepancy 
shown in Figure 13 and, in any case, are not credible because of the instability against 
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outburst (the standard model discrepancy from the roughly 6000K crossover to instability 
becomes worse and the proposed stream impact heating becomes smaller with reduced M). 
Figure 13 contrasts strongly with the very good fit of our model in Figure 10. We argue that 
no standard model synthetic spectrum, for any acceptable M, fits the observed spectrum 
and that variation of the RW Sex distance, accompanied by a compensating variation in M 
will not produce an acceptable standard model fit to the observed spectra. 

We also argue, from the above discussion of Figure 13, that variation of the WD T ef j 
cannot be used to produce an acceptable standard model fit to the observed spectrum. 
Addition of a boundary layer with a long wavelength tail in the FUV would only make the 
Figure 13 fit worse. In Figure 11 the synthetic spectrum falls below the observations near the 
FUV limit and addition of a boundary layer contribution could improve the fit. However we 
have already suggested that emission from the chromosphere could explain the discrepancy 
and we believe the evidence is too weak to support a proposal for a boundary layer. 

We are left with variation of EXP, from our previous discussion, as the only viable 
parameter to produce an acceptable fit to the observed RW Sex spectrum. 



9. The time-series FUSE spectra 



Prinja et all ( 120031 ) show that there is large orbit-to-orbit variability of the time-series 



FUSE spectra. The changes are primarily in the absorption strengths of t he high exci 



tation lines, associated with variable blue shifts tied to the orbital period. iPrinja et al. 



( 120031 ) argue that the variation arises from blueshifted absorptive changes as opposed to a 
blueshifted emission component and they suggest that the outflow wind is oblique rather 
than a symmetric bipolar wind with the oblique wind possibly seated on a warped or tilted 
disk. 

We have used our best fit model to test the tilted disk suggestion. We first subtracted 
our model synthetic spectrum, scaled to the distance of RW Sex, from the 25 individual time- 
series spectra. This step produces residuals which would display phase-wise flux variations 
if the disk is tilted. We then chose one difference spectrum (the second in the sequence) 
and successively overplotted the other difference spectra. This comparison verified the large 
case-to-case variation in the strengths of the absorption lines as well as the variable blue shift 
but there was no indication of a vertical shift in successively overplotted difference spectra. 
If the disk is warped or tilted there should be a vertical displacement of the entire difference 
spectrum that is tied to the orbital phase and we find none. Figure 14 and Figure 15 
illustrate the comparison and have the most striking differences. Superposition of the black 
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difference spectrum and other difference spectra of similar orbital phase shows close matches, 
indicating that the spectral changes primarily connect to orbital phase. 



This result agrees with the lPrinja et al.l ( 120031 ) proposal that the spectral variation arises 
from blueshifted absorptive changes. As an alternative to the the tilted disk suggestion, we 
believe a likely explanation of the non-symmetrical wind is an interaction with the debris 
from the mass transfer stream impact on the accretion disk. iLubow fc Shul ( 119761 ) find that 
the stream impact on the accretion disk leads to material flowing over th e accretion disk 
from the impact location. Hydrodynamical 2D models of the stream impact (jR6zyczkalll987l . 
19881 ) identify two shock waves: (1) a shock on a plane perpendicular to the orbital plane, 
roughly bisecting the angle between the stream and the rim and terminating at the upstream 
edge of the stream, and (2) a shock slightly mo re inclined to the stream and extending far 
into the disk. Although the simulation is 2D, iRozyczkal ( 119871 ) sta tes that a bow shock 



will d evelo p, prospectively le a ding to vertical expansion upstream. iLivio. Soker fc Dgani 



( 119861 ) and lArmitage fc Livid (119981 ) perform a 3D simulation and find that material from 
the stream flows over the disk if cooling is efficient, applicable to low M cases, and is more 
like an explosion in high M cases, leading to a bulge extending along the disk rim. 

Impact of the mass transfer stream on the accretion disk might be expected to produce 
a detectable rim hot spot, but the absence of phase- wise flux variability demonstrates that 
evidence is lacking for a sufficiently luminous rim hot spot to be detectable. This justifies 
our omission of a hot spot in the system model (§6). 



10. Discussion 



As discusssed in the Introduction, IWadd ( 119881 ) showed that neither steady state model 
accretion disks based on Planck functions nor stellar model atmospheres could simultaneously 
fit the colors and absolute luminosities of a set of NL systems (RW Sex was included in the 
study). Tomographic analyses of accretion disks ( IRutten. van Paradijs fc Tinbergenlll992[ ) 
show clear departures from standard model temperature profiles (FKR), with a spectral 
gradient typically less steep than for a st andard model. Re c ent analyses of the NL systems 
MV Lyr, IX Vel, QU Car, and UX UMa dLinnell et allbood I2007L l2008al lpn used an annulus 
model explicitly representing standard model accretion disks and showed, in each case, that a 
standard model accretion disk synthetic spectrum could not accurately fit obs erved spectra. 



On th e other hand, a standard model fits the observed spectrum in V3885 Sgr (ILinnell et al. 
20091 ). In the latter paper we speculated that a given system may at times fit a standard 



model and at other times show departures as observed in the case of UX UMa IU E spectra. 
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In the specific cases of non-standard model accretion disks listed above, the fits achieved 
depended on ad hoc changes in the temperatures of a few accretion disk annuli. The required 
departure from a standard model accretion disk temperature profile to fit a model synthetic 
spectrum to an observed spectrum has not previously been characterized analytically. This 
study shows, by detailed xl analysis, that the observed RW Sex spectrum accurately matches 
a synthetic spectrum which follows from an analyti c expre ssion, equation 3, that differs from 
the standard model. Referring back to the IWadd (119881 ) results, we conclude the problem 
is not a failure of synthetic spectra to represent the standard model but rather observed 
spectra, here RX Sex, fail to conform to the standard model. A further consequence is 
that the lower temperature gradient produces a higher rim T e g, preserving stability against 
outburst to a lower M. 

We have taken extensive precautions to consider physical effects that might otherwise 
affect the accuracy of the calculated synthetic spectra: the annulus synthetic spectra are 
non-LTE and important metals have been included as explicit ions in the calculations. It 
is curious that the empirically-determined exponent in equation 3, equal to 0.25 for the 
standard model, is closely equal to 1/2 of the standard model value for a range of M values. 
Equation 3 was not derived from a physical analysis and it is of interest to consider its 
implications. The less steep temperature profile than the standard model suggests additional 
energy deposition terms in a given ann ulus than considered in th e standard model, such as 
radial energy transfer between annuli ( IPopham &: Narayanl Il995[ ). but it seems intuitively 
doubtful that the latter effect could produce the large departure observed. In any case, 
that physical effect would be present for all accretion disks and would leave the case-to-case 
differences unexplained. Our assumption of a constant a for the entire accretion disk could 
be challenged but, as FKR discuss in connection with their equation 5.18, the energy flux 
through an accretion disk face, in the standard model, is independent of viscosity and is 
not an explicit function of a: to first order the radial temperature profile is independent 
of a. The derivation of equation 2 neglects possible magnetic effects but their inclusion, as 
mentioned in §6, currently is subject to some uncertainty. Further consideration of this topic 
is beyond the scope of this paper. 



We have neglected irradiation of the annuli by the WD. ISmak! (119891 ) shows that the 
effects are small for stationary disks. In RW Sex the presence of a disk chromosphere would 
reduce irradiation effects still further. 

It is of interest to place the results of this paper in the context of other approaches to the 



departure o 



obse rvational data from the standard model. IX Vel provides an illustration; 



Long et al.l (119941 ) consider several options including removal of all radiative flux contribution 



from inner annuli, providing a constant T e g inner region, and setting the formal inner radius 
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of the accretion disk to some multiple of the WD radius. These authors show that the 
last option, with -R rnin ~2.6-R W H , gives a reasonable fit to the HUT (Hopkins Ultraviolet 



Telescope) data. iLinnell et al.l (120071 ) showed that (for IX Vel) an accretion disk with a 
constant temperature inner section and the remainder following the standard model gives 
a good fit to FUSE and STIS spectra. Both approaches are ad hoc and the truncated 
disk model leaves unanswered the question of whether a wind can carry away all of the 
mass transfer stream, starting at the truncation radius. The ad hoc processes are localized 
perturbations of an otherwise standard model accretion disk; the analytic representation 
found here suggests existence of a process that affects the entire accretion disk. 

Our model does not include a boundary layer (BL). As discussed by FKR, half of the 
potential energy liberated in the fall of the mass transfer stream from the LI point (essen- 
tially from infinity) to the WD surface appears as radiated energy. Energy conservation 
requires an accounting for the other half and the presence of a hot BL is a common pre- 
scription. A problem is that BLs with the prescribed pro perties typically are not observed 
in high M systems or occur at unobserved wavelengths (jCordoval 119951 ). The absence of 
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In the case of the low M dwarf nova (DN) U Gem, lLong. et al.l (119961 ) used eclipse 
data during outburst to determine that an emitting region is present with a temperature of 
~140,000K, and a size approximating that of th e WD, consistent wi th a boundary layer with 
a luminosity comparable to the disk luminosity. ISzkody et al.l (119961 ) also found that in quies- 
cence U Gem has a relatively hard X-ray spectrum with the emissio n confined to a smal l area, 
supporting the interpretation of a boundary layer. More recently, iPandel et al.l ( 120051 ) show 
that, for 9 DN in quiescence, X-ray observations require a model in which a hot boundary 
layer is present and has a lumino s ity ap proximately equal to the accretion disk luminosity. 
In the cases cited by lPandel et al.l ( 20051 ) there is no boundary la yer problem and the derive d 
boundary layer structure is consistent with theoretical models (INarayan fc Pophaml 119931 ). 
The derived M rates for the DN in quiescence are M~10 -12 to 1O _11 M yr^ 1 . The M rate 
for RW Sex is between two and three dex larger. 

One proposed explanation for lack of evidence for a high luminosity boundary layer is 
that the WD is rotating clo se to Keple r ian ro tation, with no BL predicted. However, HST 
measurements starting with ISion et al.l (119941 ) have measured rotational WD velocities too 
low to explain "missing" BLs. An alternative explanatory theme is that the a bsence of an 



obser ved BL associates with a wind, and in RW Sex a wind is clearly present (IPrinja et al. 
20031 ). but the adequacy of this explanation remains undemonstrated. 
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11. Summary 

The BINSYN program suite has been used to calculate model synthetic spectra for 
RW Sex in comparison with a combined HST, FUSE and IUE spectrum. From evidence 
reported in this paper, the FU SE spectrum was obtained, apparently, at a time of reduced 
accretion disk luminosity and to achieve consistency it is necessary to divide the FUSE 
spectrum by 0.82 in combining it with the other observed spectra. 

The models include a range of M values and each model consists of 45 annuli calcu- 
lated with the Hubeny program TLUSTY; the Hubeny program SYNSPEC produces the 
synthetic spectrum for a single annulus, with a resolution of 0.1 A, includes significant met- 
als as explicit ions, and is a non-LTE model for all annuli that contribute significantly to 
the system synthetic spectrum. A Hipparcos parallax, determining D = 289pc, fixes the 
divisor that enables the system synthetic spectrum to be matched to the observed spectrum 
and establishes the fit on an absolute flux basis. Initial tests with a range of M standard 
model accretion disks show that no standard model can fit the observed spectra of RW Sex. 
New models with a generalized form of the standard model equation show that, for RW Sex, 
and based on a xt analysis, the exponent which determines the radial temperature profile is 
approximately 0.125 rather than the standard value 0.25. The optimum M for D = 289pc is 
M = 5.75xl0~ 9 M & yr" 1 , assuming the WD T cS = 50, 000K. A separate xl analysis shows 
that the data sensitivity to the WD T eff is too small to permit an actual determination of 
the WD T cS but that a value of 50,000K is a reasonable choice. 

Tests based on a separately determined D = 150pc, more recent than the Hipparcos 
determination, show that a non-standard model with M = 2.0xl0~ 9 M Q yr -1 produces a 
comparable synthetic spectrum fit to the observed spectrum as for D = 289pc (Hipparcos) 
and the xl difference in the two cases is too small to be significant. We tabulate a best fit 
model for the D = 150pc case. 

The residual spectra between individual time-series spectra (whose average is the FUSE 
spectrum described above) and the best fit synthetic spectrum model permit a study of 
the proposal for a warped or tilted accretion disk. The absence of a phase-synchronized 
vertical shift in the difference spectra rules out a tilted disk. An alternative proposal for the 
asymmetric wind is an interaction between the wind and the debris from the mass transfer 
stream impact on the accretion disk. 
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Table 1. FUV Observations of RW Sex: FUSE & HST Spectra 



Instrument 


Date 


time 


Exp. time 


Dataset 


Aperture 


Operation 


Wavelengths 




(dd/mm/yy) 


(hh:mm:ss) 


(sec) 




or Grating 


Mode 


(A) 


FUSE 


13-05-01 


14:48:05 


25,614 


B1040101 


LWRS 


TTAG 


904-1188 


HST/GHRS 


04-05-96 


13:03:53 


544 


z37v0104t 


G140L 


ACCUM 


1367-1663 


HST/GHRS 


04-05-96 


13:17:58 


544 


z37v0105t 


G140L 


ACCUM 


1367-1663 


HST/GHRS 


04-05-96 


13:30:59 


544 


z37v0106t 


G140L 


ACCUM 


1140-1435 


HST/GHRS 


04-05-96 


13:37:23 


544 


z37v0107t 


G140L 


ACCUM 


1140-1435 


HST/GHRS 


04-05-96 


14:50:24 


544 


z37v0108t 


G140L 


ACCUM 


1367-1663 


HST/GHRS 


04-05-96 


15:03:25 


435 


z37v0109t 


G140L 


ACCUM 


1367-1663 
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Table 2. FUSE ISM Line Identifications 



Ion A rcst (A) Aobs(A) Comments 



Hi 


914.29 


914.35 


18u 


Hi 


914.57 


914.58 


17u 


Hi 


914.92 


914.95 


16u 


Hi 


915.33 


915.34 


15u 


Hi 


915.82 


915.85 


14u 


Hi 


916.43 


916.47 


13u 


Oi 


916.82 


916.90 


26u 


Hi 


917.18 


917.22 


12u 


Hi 


918.13 


918.18 


llu 


Oi 


918.29 


918.48 


25u 


Oi 


918.87 


918.90 


25u 


Hi 


919.35 


919.42 


lOu 


Oi 


948.69 


948.70 


12u 


Hi 


949.74 


949.80 


4u 


Oi 


950.89 


950.90 


llu 


Pii 


961.04 


961.10 




Pn 


962.12 


962.20 




Pii 


962.57 


962.65 




Pii 


963.62 


963.65 




Pii 


963.80 


963.85 




Ni 


963.99 


964.05 


3u 


Fen 


964.30 


964.37 




Ni 


964.63 


964.70 


3u 


Pii 


964.95 


965.07 




Ni 


965.05 


965.07 


3u 


Fen 


966.20 


966.27 




Oi 


971.74 


971.80 


lOu 


Hi 


972.54 


972.65 


3u 


Oi 


976.45 


976.50 


7u 


Cm 


977.02 


977.10 


lu sharp absorption 
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Table 2 — Continued 



Ion 


Arest (A) 


Aobs(A) 


Comments 


Ol 


988.66 


988.70 


5u 


Nm 


989.80 


989.90 


lu 


Sin 


989.87 


989.90 


6u 


Si 


997.01 


997.15 




Oi 


1025.76 


1025.90 


4u 


Cn 


1036.34 


1036.35 


2u 


Cn 


1037.02 


1037.15 


2u 


Oi 


1039.23 


1039.30 


3u 


An 


1048.20 


1048.25 


2u 


An 


1066.66 


1066.70 


lu 


Nn 


1083.99 


1084.05 


lu 


Nil 


1084.57 


1084.60 


lu 


Fen 


1096.61 


1096.65 


18u 


re II 


1 nn^ to 
lUyo. (o 


lUyo.oU 


lou 


Fell 


1112.47 


1112.48 




Fell 


1121.97 


1121.9 




Fell 


1125.45 


1125.4 




Nl 


1134.17 


1134.10 


2u 


Nl 


1134.42 


1134.37 


2u 


Nl 


1134.98 


1134.92 


2u 


Fen 


1143.23 


1143.18 


lOu 


Fell 


1144.94 


1144.88 


lOu 


Pll 


1152.82 


1152.75 


3u 


? 


1181.35? 


1181.40 





Note. - In the comments column, letter u stands fo r an ultravi olet multiplet and the 
preceeding number is the multiplet number. See iMortonl (120001 120031 ). 



Table 3. FUSE Line Identifications: The blue-shifted broad absorption lines 



Ion A rcst (A) A bs(A) co-added A obs (A) diff. 



Niv 


921.46 


-920.0 


- 921.5 




921.99 


-920.0 


- 921.5 




922.52 


-920.0 


- 921.5 




923.06 


-920.0 


- 921.5 




923.22 


-920.0 


- 921.5 




924.28 


-920.0 


- 921.5 




924.91 


-920.0 


- 921.5 


Svi 


933.38 


-930. 


- 932.0 




944.52 


~940.9 


- 942.3 


Piv 


950.66 


~947.0 


- 949.2 


Cm 


977.02 


-973.5 


- 975.2 


Nm 


989.80 


-988.0 


- 988.4 


Nm 


991.51 


-988.0 


- 990.4 


Nm 


991.58 


-988.0 


- 990.4 


Ovi 


1031.91 


1028.9? 


- 1029.6 




1037.61 


1034.6? 


- 1035.9 


Siv 


1062.65 


-1060.7 


- 1061.4 


Siv 


1072.97 


-1070.4 


- 1071.8 




1073.52 


-1070.4 


- 1071.8 


Pv 


1117.98 


-1115.2 


- 1116.5 


Piv 


1118.55 


-1115.2 


- 1116.5 


Si iv 


1122.49 


-1120.0 




Pv 


1128.01 


-1125.4 


- 1126.6 


Si iv 


1128.33 


-1125.4 




Sim 


1144.31 


1142.7 






1144.96 


1142.7 






1145.11 


1142.7 






1145.18 


1142.7 






1145.67 


1142.7 




Cm 


1174.90 


1172.7 


- 1174.1 




1175.26 


1172.7 


- 1174.1 
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Table 3 — Continued 



Ion A rcst (A) A bs(A) co-added A bs(A) diff. 



1175.60 1172.7 ~ 1174.1 

1175.71 1172.7 ~ 1174.1 

1176.00 1172.7 ~ 1174.1 

1176.40 1172.7 ~ 1174.1 



-33 - 



Table 4. HST/ GHRS Line Identifications 



Ion 


. . 

Arcst(A) 


Aobs(A) 


-shift (A) 


Origin 


Qi TTT 
Ol III 


11/1/1 Q1 

1144. ol 


1 1 A o c; 
~ 114Z.0 


o c; 
-Z.o 


source 




11/1/1 CiC 

1144. yo 


i i /l o c; 
~ 114Z.0 


o c; 
-Z.O 






1 1 A C; 11 

114o.ll 


i i /l o c; 
~ 114Z.0 


o c; 
-Z.O 






1 1 /I £ 1C 

1140. lo 


i i /l o c; 
~ 114Z.0 


o c; 
-Z.o 






1 1 A £ £7 

1140.0 / 


i i /l o c; 
~ 114Z.0 


o c; 
-Z.o 




T> TT 

r II 


1 1 /1 1 so 
114Z.oy 


1 14Z.O 


n a 
-U.4 


lolVl 


13 TT 

r II 


1 1 e;o so 
lloZ.oZ 


1 1 c;o ft 
lloZ.o 


n o 
-U.Z 


lolVl 


C* TTT 

U III 


1 1 7/1 on 
iii 4.yu 


1 1 7Q /I 

~ 11 ( o.4 


O K 

-Z.o 


source 




11 f o.Zo 


1 1 7Q /I 
~ 11 i o.4 


O c; 
-Z.O 






1 1 7c; en 
11/ o.oU 


1 1 7Q A 
~ 11 1 o.4 


O c; 
-Z.O 






I 1 7c; 71 

II ( 0. / 1 


1 1 7Q /I 

~ 11 / o.4 


o c; 
-Z.O 






1 1 7fi nn 
11 / O.UU 


1 1 7Q /I 
~ 11 ( o.4 


O c; 
-Z.O 






1 1 7fi /in 
11 (0. 4U 


1 1 7Q A 
~ 11 1 o.4 


O c; 
-Z.o 




Qi TT 
Ol II 


1 1 nn /i o 
1 19U.4Z 


1 1 nn 1 
1 lyU. 1 


n q 
-U.o 


lolVl 


Qi TT 
Ol II 


1 1 no on 
llyo.Zy 


1 1 no n 
llyZ.y 


n a 
-U.4 


TQT\/T 

lolVl 


1\T T 

IN I 


1 1 nn c; c; 
llyy.oo 


i onn n 
~ 1ZUU.U 


n o 
-U.Z 


ISM, unresolved 




1 onn oo 
1ZUU.ZZ 


i onn n 
~ 1ZUU.U 


n o 
-U.Z 






1 onn 71 
1ZUU. ( 1 


i onn n 
~ 1ZUU.U 


n o 
-U.Z 




Ci TTT 
Ol III 


i onfi c; 1 
IZUb.Ol 


i one; n 
~ IZUo.U 


o n 
-Z.U 


source 




1 on7 c;o 
1ZU i .oZ 


1 one; n 
~ IZUo.U 


o n 
-Z.U 




W T 

H I 


1 01 c; ri 
IZlo.o I 


ini r O 

~ IZlo.Z 


n c; 
-U.o 


source 


N v 


1238.82 


~ 1234.2 


-4.6 


source 




1242.80 


~ 1237.9 


-4.9 




Sii 


1250.58 


1250.2 


-0.4 


ISM 


Sii 


1253.81 


1253.5 


-0.3 


ISM 


Sii 


1259.52 


1259.1 


-0.4 


ISM 


Sin 


1260.42 


1260.0 


-0.4 


ISM 


Si 


1277.22 


1277.0 


-0.2 


ISM 


Ci 


1277.26 


1277.0 


-0.3 


ISM 


Sim 


1294.55 


~ 1300.0 


-2.0 


source 




1296.73 


~ 1300.0 


-2.0 
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Table 4 — Continued 



Ion 


Arcst(A) 


Aobs(A) 


-shift (A) 


Origin 




1298.89 


~ 1300.0 


-2.0 






1298.95 


~ 1300.0 


-2.0 






1301.15 


~ 1300.0 


-2.0 






1303.32 


~ 1300.0 


-2.0 




PlI 


1301.87 


1301.7 


-0.2 


ISM 


Ol 


1302.17 


1301.7 


-0.5 


ISM 


Sin 


1304.37 


1303.9 


-0.5 


ISM 


Sim 


1328.81 


~ 1327.0 


-1.8 


source 


Cii 


1334.53 


1334.2 


-0.3 


ISM 


Cii 


1335.70 


1335.3 


-0.4 


ISM 


Si iv 


1393.76 


~ 1392.0 


-1.8 


source 




1402.77 


~ 1401.0 


-1.8 


source 


Sin 


1526.71 


1526.5 


-0.2 


ISM 


Civ 


1548.19 


~ 1542.2 


-6.0 


source 




1550.77 


~ 1544.7 


-6.1 


source 


Hen 


1640.40 


~ 1638.6 


-1.8 


source 



All the lines from the source are broad, asymmetric, blue-shifted by 2 to 6A , and belong to 
higher order ionization species, indicating that they form in a hot expanding corona above the 
disk and/or white dwarf. The lines from the ISM are sharper, blue-shifted by 0.2-0.5A , and 

belong to lower order ionization species. 
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Table 5. IUE Observations of RW Sex 



Date 


time 


Exp. time 


Description 


Aperture 


Wavelengths 


(dd/mm/yy) 


(hh:mm:ss) 


(sec) 






(A) 


30-05-78 


16:45:00 


1230 


LWR01583 


LARGE 


1851-3350 


30-11-78 


07:06:07 


714 


LWR03071 


LARGE 


1851-3350 


29-12-79 


19:44:20 


900 


LWR06494 


LARGE 


1851-3350 


30-05-78 


18:39:00 


360 


SWP01671 


LARGE 


1150-1979 


30-11-78 


08:07:12 


714 


SWP03494 


LARGE 


1150-1979 


29-12-79 


19:00:42 


900 


SWP07500 


LARGE 


1150-1979 
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Table 6. RW Sex Initial System Parameters 



parameter 


value 


parameter 


value 


M wd 


O.9OM a 


i 


34° d 


M 2 


O.674M b 


M 


LOxlO^Meyr- 1 c 


P 


0.24507 day c 


d 


289 pc f 



Beuermann. Stasiewskifc Schwopel ( 1992 ) value of 



1.35 and M 2 from iKnigg el J20061 120071 ) 
feniggel fl2006l . l2007h calibrated P : M 2 relation 



Beuermann. Stasiewskifc Schwopel (Il992l ) 



MBeuermann. Stasiewskifc Schwopel (Il992l ). i between 28 c 
and 40° 

e arbitrary initial test value 

f Hipparcos value 
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Table 7. Properties of accretion disk with mass transfer rate M = 5.0x10 9 M yr 1 and 

WD mass of 0.90M m . 



r/r wdfi 






m 


log 


g 








Ne 


Tptoss 


1 


36 


66422 


1 


21E4 


7 


07 


5 


65E7 


2 


87E17 


1 


76E4 


2 


00 


59554 


1 


48E4 


6 


80 


9 


79E7 


1 


89E17 


2 


09E4 


3 


00 


48157 


1 


44E4 


6 


50 


1 


66E8 


1 


11E17 


2 


47E4 


4 


00 


40477 


1 


35E4 


6 


29 


2 


38E8 


7 


39E16 


2 


81E4 


5 


00 


35109 


1 


25E4 


6 


12 


3 


14E8 


5 


35E16 


3 


16E4 


6 


00 


31148 


1 


17E4 


5 


98 


3 


93E8 


4 


11E16 


3 


53E4 


7 


00 


28096 


1 


10E4 


5 


86 


4 


75E8 


3 


32E16 


3 


93E4 


8 


00 


25664 


1 


04E4 


5 


75 


5 


60E8 


2 


79E16 


4 


34E4 


9 


00 


23676 


9 


83E3 


5 


66 


6 


42E8 


2 


92E16 


4 


78E4 


1U 


uu 


22016 


9 


37E3 


5 


58 


7 


28E8 


2 


55E16 


r 




2ohj4 


12 


00 


19393 


8 


60E3 


5 


44 


9 


08E8 


2 


03E16 


6 


18E4 


14 


00 


17404 


7 


98E3 


5 


32 


1 


09E9 


1 


65E16 


7 


15E4 


16 


00 


15837 


7 


47E3 


5 


21 


1 


28E9 


1 


38E16 


8 


10E4 


18 


00 


14567 


7 


04E3 


5 


12 


1 


47E9 


1 


17E16 


8 


95E4 


20 


00 


13513 


6 


67E3 


5 


03 


1 


66E9 


1 


00E16 


9 


64E4 


22 


00 


12623 


6 


35E3 


4 


96 


1 


85E9 


8 


68E15 


1 


01E5 


26 


00 


11196 


5 


82E3 


4 


82 


2 


24E9 


6 


71E15 


1 


03E5 


30 


00 


10099 


5 


39E3 


4 


70 


2 


59E9 


6 


58E15 


9 


59E4 


35 


00 


9033 


4 


97E3 


4 


57 


3 


08E9 


5 


23E15 


9 


08E4 


40 


00 


8199 


4 


62E3 


4 


45 


3 


44E9 


1 


51E14 


8 


71E4 


45 


00 


7526 


4 


33E3 


4 


24 


3 


08E9 


6 


83E13 


8 


52E4 


50 


00 


6970 


4 


09E3 


4 


07 


2 


82E9 


6 


97E13 


8 


31E4 



Note. — Each line in the t able represents a separate annulus. 
A IShakura fc Sunyaevl (119731 ) viscosity parameter a = 0.1 was 
used in calculating all annuli. The WD radius, r w d,o, is the ra- 
dius, O.OO882i? , of a zero temperature Hamada-Salpeter carbon 
model. See the text (§6.) for a discussion of the table units. 
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Table 8. Temperature profile for RW Sex accretion disk with mass transfer rate of 
M = 5.0 xKT 9 Moyr" 1 and WD mass of 0.90M Q . 



r/r wd T cff r/r wd T eff r/r wd T eff 



1 


00 


56744 


20 


71 


11994 


43 


04 


7073 


1 


14 


56744 


22 


20 


11413 


44 


52 


6900 


1 


32 


60721 


23 


69 


10893 


46 


01 


6736 


2 


85 


45281 


25 


18 


10426 


47 


50 


6582 


4 


34 


35070 


26 


66 


10004 


48 


99 


6435 


5 


83 


28942 


28 


15 


9620 


50 


48 


6296 


7 


31 


24844 


29 


64 


9269 


51 


97 


6164 


8 


80 


21891 


31 


13 


8946 


53 


46 


6038 


10 


29 


19651 


32 


62 


8649 


54 


94 


5918 


11 


78 


17885 


34 


11 


8374 


56 


43 


5804 


13 


27 


16452 


35 


59 


8119 


57 


92 


5695 


14 


76 


15264 


37 


08 


7881 


59 


41 


5590 


16 


25 


14260 


38 


57 


7659 


60 


90 


5489 


17 


73 


13399 


40 


06 


7451 


62 


39 


5393 


19 


22 


12651 


41 


55 


7256 


63 


87 


5301 



Note. - The WD radius r wd is the radius of 
a 50, 000K 0.90M M WD interpolated from Table 
4a of IPanei. Althaus fc Benvenut ol ((2000). The 
T e g = 56744K for r/r wd = 1.00 refers to the inner 
edge of the innermost annulus. The same T e g for 
J*/r w d=1.14 refers to the outer edge of the inner- 
most annulus which coincides with the inner edge 
of the next annulus. The remaining T e g values re- 
fer to the inner edge of the corresponding annuli. 
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Table 9. Values of \l 



EXP 


5.0 


5.25 


5.5 


5.75 


6.0 


6.25 


6.5 


7.0 


0.25 


86.20 


95.31 


103.75 


116.87 


130.36 


145. 


.71 


159.62 


200.61 


0.225 


79.26 


91.19 


102.89 


119.26 


136.78 


156. 


.84 


174.44 


228.35 


0.20 


63.46 


76.56 


89.58 


135.37 


128.13 


151. 


.54 


171.14 


237.15 


0.175 


39.15 


50.07 


60.87 


78.74 


98.82 


122. 


.97 


141.97 


215.97 


0.15 


19.74 


24.99 


30.39 


44.92 


61.57 


83. 


.17 


98.40 


172.34 


0.125 


23.30 


17.95 


15.51 


21.48 


30.55 


44. 


.74 


54.40 


115.13 


0.10 


55.48 


38.44 


31.36 


25.65 


25.86 


30. 


.77 


37.79 


73.20 


0.075 


95.82 


69.86 


62.42 


46.84 


40.97 


39. 


.08 


47.27 


58.82 


0.05 


124.20 


93.32 


86.65 


65.77 


57.08 


52. 


,00 


62.01 


60.47 


0.025 


137.22 


107.04 


97.42 


77.25 


67.15 


60. 


.60 


67.38 


62.87 


0.00 


140.81 


112.46 


100.65 


81.72 


71.03 


63 


.91 


66.33 


63.61 



Note. - The second and succeeding column headings are in units of 
10~ 9 M Q yr" 1 . See the text for a discussion. 
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Table 10. Values of xl 



T eS 5.0 5.25 5.5 5.75 6.0 6.25 6.5 7.0 



35 


20.51 


32.49 


40 


19.64 


17.19 


45 


21.26 


17.36 


50 


23.30 


17.95 


55 


25.03 


18.64 


60 


28.18 


20.19 


65 


30.88 


21.76 


70 


32.19 


22.62 


75 


34.53 


24.19 



22.53 18.53 15.68 

17.34 26.62 39.25 

16.22 23.85 34.71 

15.51 21.48 30.55 

15.24 19.98 27.71 

15.37 18.24 23.96 

15.90 17.38 21.61 

16.28 17.17 20.68 

17.08 16.94 19.32 



20.58 67.42 137.68 

60.51 69.78 141.41 

52.20 61.89 128.07 

44.74 54.40 115.13 

39.30 48.99 105.50 

32.49 41.43 91.71 

27.78 36.19 81.72 

25.29 33.82 76.94 

22.27 30.20 69.57 



Note. — The second and succeeding column headings are in units 
of 10~ 9 M Q jt- 1 . The T cff values are in 1000's of K. See the text for 
a discussion. 
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Table 11. RW Sex Best Fit System Parameters 



parameter 


value 


parameter 


value 


M wd 


0.90M© 


% 


34° 


EXP 


0.15 


M 2 


O.674M 


SCL 


0.51443 


M 


2.OxHT 9 M yr- 1 


T cfr (WD) 


50,000K 


P 


0.24507 day 


d 


150 pc 
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Fig. 1. — FUSE spectrum of RW Sex. The ordinate is flux in erg cm s A 
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1140 1160 1180 1200 1220 1240 1260 1280 1300 1320 




1320 1340 1360 1380 1400 1420 1440 1460 1480 




1480 1500 1520 1540 1560 1580 1600 1620 1640 1660 

Fig. 2. — HST spectrum of RW Sex. The ordinate is flux in erg cm~ 2 s^ 1 A 
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1100 1120 1140 1160 1180 1200 1220 1240 

Wavelength (A) 



Fig. 3. — Fit of FUSE, HST and IUE spectra. The black line is the FUSE spectrum as 
observed. The orange line is IUE easel; the green line is IUE case2; the cyan line is IUE 
case3. The blue line is a combination of the mean of HST z37v0104t plus z37v0105t and 
the mean of HST z37v0106t plus z37v0107t. The red line is the FUSE spectrum divided 
by 0.82. See the text for a discussion. 



E , , , i , , , , i , , , , i , , , , d 

1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 4.— Standard model synthetic spectrum fits to the observed (masked) spectrum 
of RW Sex. The standard models are for mass transfer rates of 5.0, 5.5, 6.0, 6.5, and 
7.0 xlO -9 M Q yr -1 (lowest to highest flux). See the electronic edition for a color plot. 



900 1000 1100 1200 1300 

Wavelength (A) 

Fig. 5.— As in Figure 4 but for the spectral region of the FUSE spectrum. See 
electronic edition for a color plot. 



1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 6. — Synthetic spectrum fits to observed data (black line). The orange line represents 
the value EXP = 0.25, the magenta line represents EXP = 0.125, and the red line represents 
EXP = 0.0. Note that there is little difference between the red and magenta lines at the 
longer wavelengths. 



900 1000 1100 1200 1300 

Wavelength (A) 

Fig. 7.— As in Figure 6 but for a limited range in the UV. Note the orange line (for 
EXP = 0.25) lies well above the observed spectrum (black line) and the Lyman lines are 
not deep. The red line (for EXP = 0.0) has broad Lyman lines and lies well below the 
observed spectrum. The magenta line (for EXP = 0.125) is a much better fit. 




Fig. 8. — Plot of xt values tabulated in Table 9. Standard model accretion disks correspond 
to the top horizontal line, at EXP=0.25, bounding the plot. (Compare equation 2 and 
equation 3.) Note the selection of EXP = 0.125 as roughly the best fit at all mass transfer 
rates. All models adopted a WD T eff = 50, 000K. 
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Fig. 9. — Plot of xl values tabulated in Table 10. All accretion disk models adopted a value 
EXP = 0.125. The asterisk marks the poorly constrained "best" value of T e fr, w d- As shown 
here, the originally adopted value of the WD T e g = 50, 000K was a reasonable choice. The 
lack of a closed xl contour limiting the T eff range prevents determination of the WD T eff . 
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1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 10. — Superposition of observed and model synthetic spectra. The upper black line 
is the FUSE spectrum without masking. The brown line is the IUE spectrum, without 
masking. The blue line is the mean HST z37v0104t plus z37v0105t spectrum. The cyan 
line is the mean HST z37v0106t plus z37v0107t spectrum. Note the strong chromospheric 
lines, which the model does not represent, that were masked in the solution process. The 
magenta line is the final model and the red line, almost indistinguishable from the magenta 
line but barely below it is the accretion disk contribution. The black line at the bottom is 
the 50,000K WD. 
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900 1000 1100 1200 1300 1400 1500 1600 

Wavelength (A) 

Fig. 11. — As in Figure 10 but for a restricted wavelength range. The difference between 
the magenta and red lines is barely visible. Note that this difference is the contribution of 
the 50,000K WD, shown in black at the bottom of the plot. 
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Fig. 12. — Residuals from fit of best fit model to masked observed spectrum. (Top) A 
restricted wavelength range; (Bottom) full range of combined observed spectrum. 
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1000 1500 2000 2500 3000 

Wavelength (A) 

Fig. 13. — Plot of standard model for M = 2.0 xlO -9 M Q yr -1 (gray line). The observed 
spectrum is the masked spectrum plotted in Figure 4. The contribution of the 50,000K WD 
is at the bottom. The standard model spectrum has been divided by 3.0 xlO 41 to superpose 
it on the observed spectrum. 



-55 - 



°< 



X 




900 



950 



1050 



1100 



1000 
Wavelength (A) 

Fig. 14.— Comparison of the second FUSE time-series residuals plot (black) with the 
eleventh residuals plot (magenta) at an orbital phase difference of about 0.5. Note the large 
blue shift of some of the the black plot lines and the large differences in the absorption 
line strengths. Some spectral regions match well on the two plots (950-960A, 990- 1020 A, 
1040-1055A, and 1075-llOOA), showing that there is no vertical shift of one entire difference 
spectrum relative to the other. See the text for a discussion. 
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Fig. 15. — Comparison of the second FUSE time-series residuals plot (black) with the ninth 
residuals plot (magenta) at an orbital phase difference of about 0.25. Note the large differ- 
ences in the absorption line strengths (black and magenta) and the changes from Figure 14. 
Also note that, as in Figure 14, the two difference spectra match well in the spectral intervals 
950-960A, 990-1010A, 1040-1055A, and 1075-llOOA. There is no overall vertical shift of one 
difference spectrum relative to the other. See the text for a discussion. 



